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Abstract 

In a general supersymmetric standard model there is an upper bound ruh on 
the tree level mass of the CP = +1 lightest Higgs boson which depends on 
the electroweak scale, tan (3 and the gauge and Yukawa couplings of the theory. 
When radiative corrections are included, the allowed region in the {mh,mt) 
plane depends on the scale A, below which the theory remains perturbative, 
and the supersymmetry breaking scale A^, that we fix to 1 TeV . In the minimal 
model with A = 10^^ GeV: rrih < 130 GeV and mt < 185 GeV. In non- 
minimal models with an arbitrary number of gauge singlets and A = 10^® GeV: 
rrih < 145 GeV and nit < 185 GeV. We also consider supersymmetric standard 
models with arbitrary Higgs sectors. For models whose couplings saturate the 
scale A = 10^^ GeV we find rrih < 155 GeV and rrit < 190 GeV. As one pushes 
the saturation scale A down to A^, the bounds on rrih and mt increase. For 
instance, in models with A = 10 TeV, the upper bounds for rrih and mt go to 
415 GeV and 385 GeV, respectively. 
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The most outstanding challenge for present (Tevatron, LEP) and future (LEP- 
200, NLC, LHC, SSC) colliders is the discovery of the Higgs boson which might 
confirm the standard model as the final theory of the electroweak interactions. How- 
ever, though the standard model (SM) is in excellent agreement with all precision 
measurements at present energies extensions thereof are not excluded at higher 

scales. The most appealing of these extensions, which provides a technical solution 
to the hierarchy puzzle, is the supersymmetric standard model 0. 

Supersymmetric models have well constrained Higgs sectors which can provide 
crucial tests of them. In particular, the most constraining feature of the minimal 
supersymmetric standard model (MSSM) is the existence of an absolute upper bound 
on the tree- level mass of the CP = +1 lightest Higgs boson 

rrih < mz I cos 2/3 |, (1) 

where tan/? = V2/V1, Vi = {H°). Therefore a negative result on the Higgs search 
would seem to exclude phenomenological supersymmetry at all making its search at 
future accelerators unnecessary. However relation (1) is spoiled by two effects: i) 
Radiative corrections, and ii) The enlargement of the Higgs sector in non-minimal 
supersymmetric standard models (NMSSMs). Only the simultaneous consideration 
of both effects can provide reliable bounds in general supersymmetric models. 

Radiative corrections have been computed in the MSSM by different groups using 
different methods: standard diagrammatic techniques the one- loop effective po- 
tential and the renormalization group (RG) approach including one and two 
loop corrections. All approaches provide remarkably coincident results. The latter is 
reliable provided that 

kllml, » 1, (2) 

where is the scale of supersymmetry breaking, since it amounts to a resummation 
of all leading logarithms in the effective potential. On the other hand this method is 
universal because we are assuming that the standard model holds below A^ and the 
supersymmetric theory beyond A^. Under the condition (2) supersymmetric particles 
decouple from the low-energy theory and the RG procedure is expected to provide 
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a good enough description of the radiative contribution to the hghtest Higgs boson 
mass in a general supersymmetric standard model. Our definition of lightest Higgs 
boson is the CP = +1 bosonic state whose mass is not controlled by in the sense 
that it has a finite limit when oo. In that limit supersymmetry decouples and 

the latter state becomes the SM Higgs boson. We will include radiative corrections 
for As = 1 TeV using the RG approach The radiative squared mass Ami 
/3-dependent and has to be added to the tree-level mass. 

The tree-level bound (1) does not hold in NMSSMs. The case of the MSSM plus a 
singlet with coupling \to Hi - H2 was first studied in |T0| where a tree-level bound 
was found as 



ml < \ cos' 2/5 + ^^'"^f sin^ 2/3 ) m|, (3) 



which is As-independent. From (3) we see that the bound on mh is linked to the bound 
on A if we require the theory to remain perturbative between A^ and A. For A = Aqut, 



the unification scale of gauge coupling constants, the bound (3) was studied in 111], |12 



and [0, where radiative corrections were properly included. The case of one extra 
singlet has been recently reconsidered in [Q, where the dependence of m^ on A was 
studied; |T5[, where radiative corrections where considered in the effective potential 



approach; and |T6[, where comparison with the corresponding non-supersymmetric 
scenario was established. All of these results agree, when they overlap, with our 
previous calculation |T^ within less than 5%. More general models, e.g. the MSSM 
plus any number of singlets or three SU (2) triplets (whose vacuum expectation values 
(VEVs) can respect the custodial symmetry at tree-level), were presented by ourselves 
in [jrT], |T3| and, more recently, also considered in ]T7| . 

In this paper we will present upper bounds on the lightest Higgs boson mass in 
a general class of models: supersymmetric standard models with an arbitrary Higgs 
sector. We will assume: 
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Two doublets Hi'\ H2'\ with hypercharges Y = ±1/2, coupled to quarks and 
leptons in the superpotential 

U = htQ ■ H^2^W + hQ ■ h['^D' + KL ■ Hi^'^E', (4) 

plus an arbitrary number of extra pairs Hi \ \ j = 2, ...,d + 1, decoupled 
from quarks and leptons in order to avoid dangerous flavor changing neutral 
currents [|18] . 



• Gauge singlets S^"'\ cr = 1, ...,ns- 

• SU{2) triplets S^'*), a = l,...,to, with Y = 0. 

• SU{2) triplets i = with Y = ±1. 

Notice that the above extra Higgses are the only ones that can provide renormal- 
izable couplings to h[^^ ■ H2^\ h[^^h[^^ and H^^ h!^^ in the superpotential. Other 
(more exotic) Higgs representations will only contribute to the /3-functions of the 
gauge couplings and give lower values of the upper bound nih. Since we are only 
interested in absolute upper bounds we can disregard them. Therefore the above 
class of extra Higgs fields defines the most general Higgs sector in supersymmetric 
standard models as far as the issue of putting upper bounds on their lightest Higgs 
boson mass is concerned. 

The most general renormalizable superpotential for the above Higgses can be 
written as fh + g where: 



(5) 



Summation over repeated indices is understood and traces are taken over matrix 
indices in the field decomposition: 

( I Ht \ 



(7) 
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(8) 



(9) 



(10) 



where we have dropped the multiphcity indices. 

(i) (i) 

By making a unitary transformation in j space to the Higgs doublets H^'^H^-" 
we can assume, without loss of generality, that only Hi'^ and H2'^ take a non-zero 
VEV This requires the cancellation of various Yukawa couplings in eq. (5). In 



particular we find the condition 

x'r = x'r = xi"' = xi'' = o (jvi), (11) 

that will be assumed from here on. 

After imposing condition (11) we can write fh = / + •••, where the ellipsis involves 
only fields which do not acquire VEV, i.e. (fh) = (/), as 

/ = Ai ■ SH",H° - ^ ■ CH°,m + XI ■ ViHlHl + X2 ■ ^1i^2°^2- (12) 



Where we use the notation 



(Ai)'^ = \\ 



llo- 

1 



{X2Y = >^'"' 



2 ' 



(13) 



i/° = i/p (^ = 1,2). 

Using the superpotential (12), and the fact that the smallest eigenvalue of a real, 
symmetric nxn matrix is smaller than the smallest eigenvalue of the upper left 2x2 
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submatrix, one can easily prove that the hghtest Higgs boson mass has an upper 
bound given by |Tl[ 



ml/v^ <li9' + g") cos2 2/5 + {t +^ a') sin^ 2/3+ xl cos^ /5+ xl sin^ /3, (14) 

where v"^ = vf + v\ and g,g' are the SU{2) x f/(l)y gauge couphngs. In particular 
the bound (3) is recovered when A2 = Xi = X2 = and the bound (1) in the MSSM 
when also Ai = 0. 



As it has been repeatedly noticed [0, [T^ the bound (14) is independent of the soft- 



supersymmetry breaking parameters, or any supersymmetric mass terms. It is only 
controlled by v and dimensionless parameters ((7, (7', tan /5, and Yukawa couplings). 
Since the former is fixed by the electroweak scale, the latter will determine the bound 
(14). In particular the upper bound on the right hand side of (14) comes from 
the requirement that the supersymmetric theory remains perturbative below A. To 
guarantee this condition we need to solve the renormalization group equations (RGEs) 
of all gauge and Yukawa couplings of the theory. 

In the general supersymmetric theory defined by the superpotential (4-6) with the 
condition (11), the relevant RGEs to the bound (14) can be written, for the Yukawa 
couplings as: 



4vr2| A? = [-Ig' - \g'^ + f (/^? + hi) + ^ % +3(x? + xl)+2 \{) x{ 

+ E (A^^A^')^ + ^'^^ 

k,l^l ^ IJ-,f 



^ k,l^l ^ k,l ^ k,l 

47r2| xl = {-h' - y + ^hl + l t +7 xl + A?) xl 
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2^^« _ ( 3^2 13^/2 8^2 I Q/,2 I 1 !,2 I 3 I Q -2 1 r2 



^ = -1^?^ - 11^7'^ - %gl + 3/.,^ + i/^^ + I \2 +3 X2 +i Ai U 



it 



+ 1 \ Y.>^r\r + ^Y.xrxr + ^ E^f ^Af ^^^^ 



2(i/?-6 _ / 3 2 7 ^/2 8^2 I qa2 1 1 a2 1 3 I Q 1 1 



^'-^ = - 1^^?'^ - 1^?.^ + + i/i,^ + I A2 +3 x\ +i Ai U 



+ f z E^2'"^2''^ + 3Exf ^x}'^^ + ^ Y.>^rxr\ h, 



(20) 



and for the gauge couplings as: 

16n^^ ={l + 2to + 4h + d)g^ (21) 

16vr2^ =(U + 6t, + d)g'^ (22) 

167r2% =-3^7f, (23) 

where gg is the SU{3) gauge couphng Notice that condition (11) is stable under 
the RGEs. 

The r- Yukawa coupling, hr, will be neglected as compared to hb, since hr/hb = 
mT-/mf,. The bottom Yukawa coupling can be important for tan/? ^ 1 and will be 
kept along with the top Yukawa coupling. They are given by 

(24) 

V2 mw 



^For simplicity we did not present the explicit RGEs for the couplings Af''^, X^J"^, Xi^ , xi^ ^ ^ajk, 

f^abc cind Ka-fii" 
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and fixed by the boundary conditions: mt{2mt) = rrit and mb{2mb) = 5 GeV. For 
the gauge couplings we will take the boundary conditions: 

aEMiMz) = j^, smHwiMz) = 0.23, a,(M^) = 0.12 . (25) 

The Yukawa couplings involved in (14) will be let to acquire any perturbative 
value maximizing the bound (14). In general the bound (14) is maximized whenever 
some of the involved couplings saturate the scale A. A particular coupling A (gauge 
or Yukawa) is said to saturate a scale A if 

A'(Q')/47r < 1, (26) 

for < A^, and the equality in (26) holds for = A^. 

Once we have the set of RGBs (15-23) we can systematically analyze the case 
of different supersymmetric standard models characterized by different Higgs sec- 
tors. The different cases are parametrized by the number of Higgs representations 
{ns,d,to,ti). For simplicity we will assume to = ti = t such that the custodial sym- 
metry (and so the tree level value p = 1) can be respected by means of a very simple 
relation between the VEVs of S, \E'i and \l/2. In this way the gauge couplings will be 
parametrized by the single parameter 

N = 6t + d, (27) 

and the bound will depend on Ug and A^. In all cases the bound will be a function of 
A, and so it is important to discuss the criteria to fix it. 

• In theories with gauge coupling unification at the scale Aqut, the natural choice 
is fixing A = Aqut- The case satisfying this condition is the MSSM plus an 
arbitrary number of gauge singlets, i.e. N = 0, arbitrary. Of course putting 

= we recover the MSSM. 

• In theories without gauge coupling unification, 7^ 0, we will assign to any 

scale A the value of A^ such that the gauge coupling g saturates it 0. Once we 
^We can see from (21-23) and conditions (25-27) that every scale A is always saturated by g: this 
means that g is the first gauge coupling to go non-perturbative. 
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have fixed A^, the bound (14) is maximized for Ug ^ 0, and so we will assume 
there are gauge singlets. The actual value of the bound does not depend on the 
particular value of Ug, provided that > 0. 

The key observation to maximize the bound (14) is to notice that Aj and Xi (^ = 
1,2) are maximized by the values 

Ar = Ar = xf = xf = 

(28) 

^ajk l^abc ^ajxv 0, 

since they contribute to the renormalization of and xi {i — 1,2) with positive 
definite terms. Eq. (28) is stable under the RGEs, so we can impose it in order to 
obtain upper bounds on rrih. It might happen in specific models that we would need 
to departure from condition (28) to encompass the low energy experimental bounds. 
For instance, this could be the case when some of the couplings involved in (28) are 
necessary to explicitly break a global symmetry; in that case, to avoid a massless 
axion, we should put them to non-zero values. However, since we are dealing only 
with absolute upper bounds we will use condition (28) in the rest of this paper. 

We have analyzed different cases characterized by different values of A. In Fig. 1 
we plot nih as a function of for the model N — 0, Ug > (i.e. the MSSM enlarged 
with Higgs singlets) and different values of tan /3. In Fig. 2 we plot the case saturating 
the scale A = 10^^ GeV, which corresponds to A^ = 5. The model with t — l,d = 
{N = 6), which saturates the scale A = 10^^ GeV, is shown in Fig. 3. We see that the 
value of rrih increases as we let the saturation scale A to go down to A^. In particular 
we show in Fig. 4 the plot corresponding to A = 10^° GeV {N — 10). Finally we show 
in Fig. 5 the allowed region in the {mh,mt) plane for the MSSM (long-dashed); the 
model N = 0,ns > with A = Aqut (short-dashed), and different models saturating 
different scales, from A = 10^^ GeV to A = 10 TeV (solid). 

In conclusion, we have computed upper bounds on the mass of the lightest Higgs 
boson in a general class of supersymmetric standard models, characterized by arbi- 
trary Higgs sectors. For a given scale A the corresponding upper bound is saturated 



8 



by the model whose gauge couphngs become non-perturbative at that scale. We have 
taken the scale of supersymmetry breaking = 1 TeV and assumed a supersymmet- 
ric theory between A and A^ and the standard model below A^. In detailed studies 
of particular models, threshold effects of supersymmetric particles should be taken 
into account, as well as physical conditions on the Yukawa couplings that we are not 
considering {e.g. unbroken color and/or electric charge). These effects could slightly 
change our absolute upper bounds. Radiative corrections are taken into account us- 
ing the (universal) renormalization group method. Diagrammatic techniques should 
give similar results provided that A^ ^ tti"^. (This has been explicitly checked in the 
simplest non-minimal model with a singlet.) For lower values of As the detailed calcu- 
lation should be done for every particular model. We did not analyze the dependence 
of the upper bounds on A^ but we expect it not to be too dramatic provided A^ is 
kept inside its phenomenological range. For all the above reasons our results should 
be taken as rough estimates to guide the eye in supersymmetric standard models. 
They are strongly dependent on the scale A. In fact for A from 10^^ GeV to 10'^ GeV 
the upper bounds on range from 130 to 415 GeV, and the bounds on nit from 185 
to 385 GeV. For A = A^ the bound coincides with that in the non-supersymmetric 
standard model |pO|| . 

We have not considered in this paper other possible generalizations of the MSSM 
that could, in turn, modify the upper bounds on the lightest Higgs boson mass. One 
of them is the introduction of extra colored multiplets, as e.g. extra (fractions of) 
generations. They would provide large radiative corrections, through the Yukawa 
couplings of the additional up quarks, to the Higgs boson mass for the case of a not- 
so-heavy top quark Another possible generalization is the presence of an extra 
gauge group G, with couplings ga and generators which would amount to a 

correction to the tree-level bound (14) as 

Ami = E (^f cos^ (3 - sin' (3) ' (29) 

a 

where T^^ = (H^^^T'^hI^^) /vf. In that case the Yukawa couplings in (4) and (5) should 
respect the gauge symmetry G and all the bounds could change dramatically. 
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Figure Captions 

Fig. 1 Upper bounds on the lightest scalar Higgs boson in NMSSM with singlets 
and different values of tan /3 (sohd) . 

Fig. 2 The same as in Fig.l but for a model saturating A = 10^^ GeV. 

Fig. 3 The same as in Fig.l but for a model saturating A = 10^^ GeV. 

Fig. 4 The same as in Fig.l but for a model saturating A = 10^^ GeV. 

Fig. 5 Allowed regions in the (m/i,mt) plane for different supersymmetric stan- 
dard models. The solid curves correspond to models saturating the scales 
A = 10^^ 10l^ 10^°, 10^ 10^ and 10^ GeV. 
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